Chromosome-level Genome Assembly of the Yeast Lodderomyces beijingensis Reveals the Genetic Nature of Metabolic Adaptations and Identifies Subtelomeres as Hotspots for Amplification of Mating Type Loci.

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterize its genome sequence.

y-mtPTM: Yeast mitochondrial posttranslational modification database.

One powerful strategy of how to increase the complexity of cellular proteomes is through posttranslational modifications (PTMs) of proteins. Currently, there are ∼400 types of PTMs, the different combinations of which yield a large variety of protein isoforms with distinct biochemical properties. Although mitochondrial proteins undergoing PTMs were identified nearly 6 decades ago, studies on the roles and extent of PTMs on mitochondrial functions lagged behind the other cellular compartments. The application of mass spectrometry for the characterization of the mitochondrial proteome as well as for the detection of various PTMs resulted in the identification of thousands of amino acid positions that can be modified by different chemical groups. However, the data on mitochondrial PTMs are scattered in several data sets, and the available databases do not contain a complete list of modified residues. To integrate information on PTMs of the mitochondrial proteome of the yeast Saccharomyces cerevisiae, we built the yeast mitochondrial posttranslational modification (y-mtPTM) database (http://compbio.fmph.uniba.sk/y-mtptm/). It lists nearly 20,000 positions on mitochondrial proteins affected by ∼20 various PTMs, with phosphorylated, succinylated, acetylated, and ubiquitylated sites being the most abundant. A simple search of a protein of interest reveals the modified amino acid residues, their position within the primary sequence as well as on its 3D structure, and links to the source reference(s). The database will serve yeast mitochondrial researchers as a comprehensive platform to investigate the functional significance of the PTMs of mitochondrial proteins.

Chromosome-Level Genome Assembly of the Yeast Candida verbasci.

Candida verbasci is an anamorphic ascomycetous yeast. We report the genome sequence of its type strain, 11-1055 (CBS 12699). The nuclear genome assembly consists of seven chromosome-sized contigs with a total size of 12.1 Mbp and has a relatively low G+C content (28.1%).

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex.

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.

Transcriptome and proteome profiling reveals complex adaptations of Candida parapsilosis cells assimilating hydroxyaromatic carbon sources.

Many fungal species utilize hydroxyderivatives of benzene and benzoic acid as carbon sources. The yeast Candida parapsilosis metabolizes these compounds via the 3-oxoadipate and gentisate pathways, whose components are encoded by two metabolic gene clusters. In this study, we determine the chromosome level assembly of the C. parapsilosis strain CLIB214 and use it for transcriptomic and proteomic investigation of cells cultivated on hydroxyaromatic substrates. We demonstrate that the genes coding for enzymes and plasma membrane transporters involved in the 3-oxoadipate and gentisate pathways are highly upregulated and their expression is controlled in a substrate-specific manner. However, regulatory proteins involved in this process are not known. Using the knockout mutants, we show that putative transcriptional factors encoded by the genes OTF1 and GTF1 located within these gene clusters function as transcriptional activators of the 3-oxoadipate and gentisate pathway, respectively. We also show that the activation of both pathways is accompanied by upregulation of genes for the enzymes involved in ß-oxidation of fatty acids, glyoxylate cycle, amino acid metabolism, and peroxisome biogenesis. Transcriptome and proteome profiles of the cells grown on 4-hydroxybenzoate and 3-hydroxybenzoate, which are metabolized via the 3-oxoadipate and gentisate pathway, respectively, reflect their different connection to central metabolism. Yet we find that the expression profiles differ also in the cells assimilating 4-hydroxybenzoate and hydroquinone, which are both metabolized in the same pathway. This finding is consistent with the phenotype of the Otf1p-lacking mutant, which exhibits impaired growth on hydroxybenzoates, but still utilizes hydroxybenzenes, thus indicating that additional, yet unidentified transcription factor could be involved in the 3-oxoadipate pathway regulation. Moreover, we propose that bicarbonate ions resulting from decarboxylation of hydroxybenzoates also contribute to differences in the cell responses to hydroxybenzoates and hydroxybenzenes. Finally, our phylogenetic analysis highlights evolutionary paths leading to metabolic adaptations of yeast cells assimilating hydroxyaromatic substrates.

Monoclonal antibodies targeting two immunodominant epitopes on the Spike protein neutralize emerging SARS-CoV-2 variants of concern.

BACKGROUND: The emergence of new SARS-CoV-2 variants of concern B.1.1.7 (Alpha), B.1.351 (Beta), P.1 (Gamma) and B.1.617.2 (Delta) that harbor mutations in the viral S protein raised concern about activity of current vaccines and therapeutic antibodies. Independent studies have shown that mutant variants are partially or completely resistant against some of the therapeutic antibodies authorized for emergency use. METHODS: We employed hybridoma technology, ELISA-based and cell-based S-ACE2 interaction assays combined with authentic virus neutralization assays to develop second-generation antibodies, which were specifically selected for their ability to neutralize the new variants of SARS-CoV-2. FINDINGS: AX290 and AX677, two monoclonal antibodies with non-overlapping epitopes, exhibit subnanomolar or nanomolar affinities to the receptor binding domain of the viral Spike protein carrying amino acid substitutions N501Y, N439K, E484K, K417N, and a combination N501Y/E484K/K417N found in the circulating virus variants. The antibodies showed excellent neutralization of an authentic SARS-CoV-2 virus representing strains circulating in Europe in spring 2020 and also the variants of concern B.1.1.7 (Alpha), B.1.351 (Beta) and B.1.617.2 (Delta). In addition, AX677 is able to bind Omicron Spike protein just like the wild type Spike. The combination of the two antibodies prevented the appearance of escape mutations of the authentic SARS-CoV-2 virus. Prophylactic administration of AX290 and AX677, either individually or in combination, effectively reduced viral burden and inflammation in the lungs, and prevented disease in a mouse model of SARS-CoV-2 infection. INTERPRETATION: The virus-neutralizing properties were fully reproduced in chimeric mouse-human versions of the antibodies, which may represent a promising tool for COVID-19 therapy. FUNDING: The study was funded by AXON Neuroscience SE and AXON COVIDAX a.s.

Surveillance of SARS-CoV-2 lineage B.1.1.7 in Slovakia using a novel, multiplexed RT-qPCR assay.

The emergence of a novel SARS-CoV-2 B.1.1.7 variant sparked global alarm due to increased transmissibility, mortality, and uncertainty about vaccine efficacy, thus accelerating efforts to detect and track the variant. Current approaches to detect B.1.1.7 include sequencing and RT-qPCR tests containing a target assay that fails or results in reduced sensitivity towards the B.1.1.7 variant. Since many countries lack genomic surveillance programs and failed assays detect unrelated variants containing similar mutations as B.1.1.7, we used allele-specific PCR, and judicious placement of LNA-modified nucleotides to develop an RT-qPCR test that accurately and rapidly differentiates B.1.1.7 from other SARS-CoV-2 variants. We validated the test on 106 clinical samples with lineage status confirmed by sequencing and conducted a country-wide surveillance study of B.1.1.7 prevalence in Slovakia. Our multiplexed RT-qPCR test showed 97% clinical sensitivity and retesting 6,886 SARS-CoV-2 positive samples obtained during three campaigns performed within one month, revealed pervasive spread of B.1.1.7 with an average prevalence of 82%. Labs can easily implement this test to rapidly scale B.1.1.7 surveillance efforts and it is particularly useful in countries with high prevalence of variants possessing only the ΔH69/ΔV70 deletion because current strategies using target failure assays incorrectly identify these as putative B.1.1.7 variants.

Nanopore sequencing of SARS-CoV-2: Comparison of short and long PCR-tiling amplicon protocols.

Surveillance of the SARS-CoV-2 variants including the quickly spreading mutants by rapid and near real-time sequencing of the viral genome provides an important tool for effective health policy decision making in the ongoing COVID-19 pandemic. Here we evaluated PCR-tiling of short (~400-bp) and long (~2 and ~2.5-kb) amplicons combined with nanopore sequencing on a MinION device for analysis of the SARS-CoV-2 genome sequences. Analysis of several sequencing runs demonstrated that using the long amplicon schemes outperforms the original protocol based on the 400-bp amplicons. It also illustrated common artefacts and problems associated with PCR-tiling approach, such as uneven genome coverage, variable fraction of discarded sequencing reads, including human and bacterial contamination, as well as the presence of reads derived from the viral sub-genomic RNAs.

The yeast mitochondrial succinylome: Implications for regulation of mitochondrial nucleoids.

Acylation modifications, such as the succinylation of lysine, are post-translational modifications and a powerful means of regulating protein activity. Some acylations occur nonenzymatically, driven by an increase in the concentration of acyl group donors. Lysine succinylation has a profound effect on the corresponding site within the protein, as it dramatically changes the charge of the residue. In eukaryotes, it predominantly affects mitochondrial proteins because the donor of succinate, succinyl-CoA, is primarily generated in the tricarboxylic acid cycle. Although numerous succinylated mitochondrial proteins have been identified in Saccharomyces cerevisiae, a more detailed characterization of the yeast mitochondrial succinylome is still lacking. Here, we performed a proteomic MS analysis of purified yeast mitochondria and detected 314 succinylated mitochondrial proteins with 1763 novel succinylation sites. The mitochondrial nucleoid, a complex of mitochondrial DNA and mitochondrial proteins, is one of the structures whose protein components are affected by succinylation. We found that Abf2p, the principal component of mitochondrial nucleoids responsible for compacting mitochondrial DNA in S. cerevisiae, can be succinylated in vivo on at least thirteen lysine residues. Abf2p succinylation in vitro inhibits its DNA-binding activity and reduces its sensitivity to digestion by the ATP-dependent ScLon protease. We conclude that changes in the metabolic state of a cell resulting in an increase in the concentration of tricarboxylic acid intermediates may affect mitochondrial functions.

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020.

SARS-CoV-2 mutants carrying the ∆H69/∆V70 deletion in the amino-terminal domain of the Spike protein emerged independently in at least six lineages of the virus (namely, B.1.1.7, B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375). We analyzed SARS-CoV-2 samples collected from various regions of Slovakia between November and December 2020 that were presumed to contain B.1.1.7 variant due to drop-out of the Spike gene target in an RT-qPCR test caused by this deletion. Sequencing of these samples revealed that although in some cases the samples were indeed confirmed as B.1.1.7, a substantial fraction of samples contained another ∆H69/∆V70 carrying mutant belonging to the lineage B.1.258, which has been circulating in Central Europe since August 2020, long before the import of B.1.1.7. Phylogenetic analysis shows that the early sublineage of B.1.258 acquired the N439K substitution in the receptor-binding domain (RBD) of the Spike protein and, later on, also the deletion ∆H69/∆V70 in the Spike N-terminal domain (NTD). This variant was particularly common in several European countries including the Czech Republic and Slovakia but has been quickly replaced by B.1.1.7 early in 2021.

Genome analysis of Candida subhashii reveals its hybrid nature and dual mitochondrial genome conformations.

Candida subhashii belongs to the CUG-Ser clade, a group of phylogenetically closely related yeast species that includes some human opportunistic pathogens, such as Candida albicans. Despite being present in the environment, C. subhashii was initially described as the causative agent of a case of peritonitis. Considering the relevance of whole-genome sequencing and analysis for our understanding of genome evolution and pathogenicity, we sequenced, assembled and annotated the genome of C. subhashii type strain. Our results show that C. subhashii presents a highly heterozygous genome and other signatures that point to a hybrid ancestry. The presence of functional pathways for assimilation of hydroxyaromatic compounds goes in line with the affiliation of this yeast with soil microbial communities involved in lignin decomposition. Furthermore, we observed that different clones of this strain may present circular or linear mitochondrial DNA. Re-sequencing and comparison of strains with differential mitochondrial genome topology revealed five candidate genes potentially associated with this conformational change: MSK1, SSZ1, ALG5, MRPL9 and OYE32.

OCT1 - a yeast mitochondrial thiolase involved in the 3-oxoadipate pathway.

The 3-oxoacyl-CoA thiolases catalyze the last step of the fatty acid ß-oxidation pathway. In yeasts and plants, this pathway takes place exclusively in peroxisomes, whereas in animals it occurs in both peroxisomes and mitochondria. In contrast to baker's yeast Saccharomyces cerevisiae, yeast species from the Debaryomycetaceae family also encode a thiolase with predicted mitochondrial localization. These yeasts are able to utilize a range of hydroxyaromatic compounds via the 3-oxoadipate pathway the last step of which is catalyzed by 3-oxoadipyl-CoA thiolase and presumably occurs in mitochondria. In this work, we studied Oct1p, an ortholog of this enzyme from Candida parapsilosis. We found that the cells grown on a 3-oxoadipate pathway substrate exhibit increased levels of the OCT1 mRNA. Deletion of both OCT1 alleles impairs the growth of C. parapsilosis cells on 3-oxoadipate pathway substrates and this defect can be rescued by expression of the OCT1 gene from a plasmid vector. Subcellular localization experiments and LC-MS/MS analysis of enriched organellar fraction-proteins confirmed the presence of Oct1p in mitochondria. Phylogenetic profiling of Oct1p revealed an intricate evolutionary pattern indicating multiple horizontal gene transfers among different fungal groups.

Shared evolutionary footprints suggest mitochondrial oxidative damage underlies multiple complex I losses in fungi.

Oxidative phosphorylation is among the most conserved mitochondrial pathways. However, one of the cornerstones of this pathway, the multi-protein complex NADH : ubiquinone oxidoreductase (complex I) has been lost multiple independent times in diverse eukaryotic lineages. The causes and consequences of these convergent losses remain poorly understood. Here, we used a comparative genomics approach to reconstruct evolutionary paths leading to complex I loss and infer possible evolutionary scenarios. By mining available mitochondrial and nuclear genomes, we identified eight independent events of mitochondrial complex I loss across eukaryotes, of which six occurred in fungal lineages. We focused on three recent loss events that affect closely related fungal species, and inferred genomic changes convergently associated with complex I loss. Based on these results, we predict novel complex I functional partners and relate the loss of complex I with the presence of increased mitochondrial antioxidants, higher fermentative capabilities, duplications of alternative dehydrogenases, loss of alternative oxidases and adaptation to antifungal compounds. To explain these findings, we hypothesize that a combination of previously acquired compensatory mechanisms and exposure to environmental triggers of oxidative stress (such as hypoxia and/or toxic chemicals) induced complex I loss in fungi.

Step-by-Step Evolution of Telomeres: Lessons from Yeasts.

In virtually every eukaryotic species, the ends of nuclear chromosomes are protected by telomeres, nucleoprotein structures counteracting the end-replication problem and suppressing recombination and undue DNA repair. Although in most cases, the primary structure of telomeric DNA is conserved, there are several exceptions to this rule. One is represented by the telomeric repeats of ascomycetous yeasts, which encompass a great variety of sequences, whose evolutionary origin has been puzzling for several decades. At present, the key questions concerning the driving force behind their rapid evolution and the means of co-evolution of telomeric repeats and telomere-binding proteins remain largely unanswered. Previously published studies addressed mostly the general concepts of the evolutionary origin of telomeres, key properties of telomeric proteins as well as the molecular mechanisms of telomere maintenance; however, the evolutionary process itself has not been analyzed thoroughly. Here, we aimed to inspect the evolution of telomeres in ascomycetous yeasts from the subphyla Saccharomycotina and Taphrinomycotina, with special focus on the evolutionary origin of species-specific telomeric repeats. We analyzed the sequences of telomeric repeats from 204 yeast species classified into 20 families and as a result, we propose a step-by-step model, which integrates the diversity of telomeric repeats, telomerase RNAs, telomere-binding protein complexes and explains a propensity of certain species to generate the repeat heterogeneity within a single telomeric array.

Mitochondrial protein phosphorylation in yeast revisited.

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

Design and Synthesis of Pyrano[3,2-b]indolones Showing Antimycobacterial Activity.

Latent Mycobacterium tuberculosis infection presents one of the largest challenges for tuberculosis control and novel antimycobacterial drug development. A series of pyrano[3,2-b]indolone-based compounds was designed and synthesized via an original eight-step scheme. The synthesized compounds were evaluated for their in vitro activity against M. tuberculosis strains H37Rv and streptomycin-starved 18b (SS18b), representing models for replicating and nonreplicating mycobacteria, respectively. Compound 10a exhibited good activity with MIC99 values of 0.3 and 0.4 µg/mL against H37Rv and SS18b, respectively, as well as low toxicity, acceptable intracellular activity, and satisfactory metabolic stability and was selected as the lead compound for further studies. An analysis of 10a-resistant M. bovis mutants disclosed a cross-resistance with pretomanid and altered relative amounts of different forms of cofactor F420 in these strains. Complementation experiments showed that F420-dependent glucose-6-phosphate dehydrogenase and the synthesis of mature F420 were important for 10a activity. Overall these studies revealed 10a to be a prodrug that is activated by an unknown F420-dependent enzyme in mycobacteria.

Mitochondrial HMG-Box Containing Proteins: From Biochemical Properties to the Roles in Human Diseases.

Mitochondrial DNA (mtDNA) molecules are packaged into compact nucleo-protein structures called mitochondrial nucleoids (mt-nucleoids). Their compaction is mediated in part by high-mobility group (HMG)-box containing proteins (mtHMG proteins), whose additional roles include the protection of mtDNA against damage, the regulation of gene expression and the segregation of mtDNA into daughter organelles. The molecular mechanisms underlying these functions have been identified through extensive biochemical, genetic, and structural studies, particularly on yeast (Abf2) and mammalian mitochondrial transcription factor A (TFAM) mtHMG proteins. The aim of this paper is to provide a comprehensive overview of the biochemical properties of mtHMG proteins, the structural basis of their interaction with DNA, their roles in various mtDNA transactions, and the evolutionary trajectories leading to their rapid diversification. We also describe how defects in the maintenance of mtDNA in cells with dysfunctional mtHMG proteins lead to different pathologies at the cellular and organismal level.

Role of folding kinetics of secondary structures in telomeric G-overhangs in the regulation of telomere maintenance in Saccharomyces cerevisiae.

The ends of eukaryotic chromosomes typically contain a 3' ssDNA G-rich protrusion (G-overhang). This overhang must be protected against detrimental activities of nucleases and of the DNA damage response machinery and participates in the regulation of telomerase, a ribonucleoprotein complex that maintains telomere integrity. These functions are mediated by DNA-binding proteins, such as Cdc13 in Saccharomyces cerevisiae, and the propensity of G-rich sequences to form various non-B DNA structures. Using CD and NMR spectroscopies, we show here that G-overhangs of S. cerevisiae form distinct Hoogsteen pairing-based secondary structures, depending on their length. Whereas short telomeric oligonucleotides form a G-hairpin, their longer counterparts form parallel and/or antiparallel G-quadruplexes (G4s). Regardless of their topologies, non-B DNA structures exhibited impaired binding to Cdc13 in vitro as demonstrated by electrophoretic mobility shift assays. Importantly, whereas G4 structures formed relatively quickly, G-hairpins folded extremely slowly, indicating that short G-overhangs, which are typical for most of the cell cycle, are present predominantly as single-stranded oligonucleotides and are suitable substrates for Cdc13. Using ChIP, we show that the occurrence of G4 structures peaks at the late S phase, thus correlating with the accumulation of long G-overhangs. We present a model of how time- and length-dependent formation of non-B DNA structures at chromosomal termini participates in telomere maintenance.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends.

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.